Beyond the time independent mean field theory for nuclear and atomic reactions: 
Inclusion of particle-hole correlations in a generalized random phase approximation 
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The time independent mean field metliod (TIMF) for scattering defines biortfionormal sets of single- 
particle wave functions and corresponding creation and annihilation operators. 2p — 2h correlations 
can be introduced through a generalized random phase approximation; Ip— Ih contributions vanish 
(Brillouin theorem). While the general variational method for scattering by Giraud and Nagarajan 
solves inhomogeneous Euler equations by inversion of the standard, hermitean Hamiltonian, the 
present approach diagonalizes a non-hermitean Hamiltonian, which carries the information about 
entrance and exit channels. 



The time independent mean field theory of collisions 
(TIMF) [Q has been tested successfully for solvable cases 
of 3 and 4 particles and has been applied to various 
nuclear and atomic reactions . It is then natural to go 
beyond the mean field approach by introducing particle- 
hole (ph) correlations in analogy to the case of bound 
states. In this letter we outline how this goal can be 
reached by a generalized random phase approximation. 
Our approach can be viewed as an alternative of more 
traditional inclusions of correlations into scattering theo- 
ries, such as the antisymmctrization of cluster structures 
[4| and fiuctuations into time dependent Hartree-Fock 
[5|, for instance. 

We start from a time independent functional Q which 
calculates the Green function D{z) between incoming 
and outgoing channel waves x,x'- The variation of 

F(*', *) x'l* >< > / < *'|(z - i/)|* >, (1) 

with Im z > gives, in appropriate normalization, 

{z - H)\^ >^ \X > , <^'\iz-H)=<x'\, (2) 

Diz) =< x'\{z - H)-'\x >-< x'l* >=< ^'\x > , (3) 

with ^',^1^' as stationary (saddle!) points of (1). For 
simplicity we assume that x ^^id x' *re Slater determi- 
nants made of N square integrable orbitals Xi ^-iid Xi: 
respectively. H is the standard full Hamiltonian of N 
particles. In the same way one may obtain the T-matrix 
by multiplying Xj x' by prior and post potentials V, V' 
in the above equations. While one normally calculates 
D as function of (complex) energy z, we now reverse the 
strategy: We rewrite (2) with (3) as 



{z-H)\^> =A|x><x'l*>, 
< ^-'Kz - i?) = A < ^-'Ix >< x'l , 



which we interpret as right and left eigenvalue equations 
of the non-hermitean Hamiltonian 



H' = H + A\x><x'\; A = 1/L> complex. 



(5) 



with eigenvalues z = z{A). Choosing fixed values of A 
in (4) and (5) preserves the equivalence of (2) and (4), 
since Eqs.(4) are invariant under normalization of 5*, ^P' 
and their solutions are stationary points of the functional 
(1) which also is invariant under normalization of ^I''. 

The discrete eigenvalues of H' are solutions of the ob- 
vious quantization condition 



=<x'\iz-H)-'\x>. 



(6) 



which holds for both right and left eigenstates according 
to Eqs.(4). It is clear that any real E of the discrete 
spectrum of H generates a complex z close to E, at least 
whenever A is small. Conversely, there could be discrete 
eigenvalues z of H' which are not in the vicinity of dis- 
crete eigenvalues oi H. As long as A is the reciprocal of 
a physical amplitude D, an additional, discrete, complex 
eigenvalue emerges, close to the physical scattering en- 
ergy E. The continuum of H' is the same as that of H 
because of the compactness of |x >< x'l- 

Square integrable left and right eigenstates of H' can 
be found from the Rayleigh-Ritz-like functional 



F'{^', *) =< *'|(z - H')\^ > , 



(7) 



(4a) 
(4b) 



where z appears as a norm and phase control Lagrange 
multiplier, actually equal to the desired eigenvalue of H', 
and the trial functions 5*, '5' are square integrable as long 
as Imz 0. Indeed, the variational equations of (7) 
induce the right and left eigenvalue equations (4). It 
should be stressed here that, from the point of view of 
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the variational principle for F' , one looks for z when A 
is given, while from the point of view of F, one looks for 
A when z is given. 

In this paper we investigate solutions of (4) where 
Re 2; can be interpreted as an energy of the continuum 
of H, with Imz > 0, and we consider 2; as a function 
of A. In other words, we seek an energy z for which 
the Green function D{z) is equal to a given value A~^. 
Our way of solving (6) for z consists in diagonalizing 
H' via suitable approximations of vl>, ^t' in the functional 
(7). This involves first a Hartree-Fock approximation 
which induces a zero order answer zo{A.), then possibly a 
Tamm-Dancoff approach (TDA) to the spectrum of H' 
and finally random phase (RPA), second order correc- 
tions Az{A), with the final result 2:2 (A) = zo{A) + Az{A}. 
First order corrections are automatically cancelled by the 
Hartree-Fock method (generalized Brillouin theorem). 

The Hartree-Fock approximation to the diagonaliza- 
tion of H' consists in restricting ^, ^' to just Slater de- 
terminants made of orbitals (pi, = 1,2,... A''. 
To simplify the corresponding variational equations, we 
redefine Xi, Xi by the conditions 



< f'i\X3 >= =< XilVj > tor i^j. 



(8) 



This just means a linear rearrangement of the orbitals 

Xi^Xi^ leaving the determinants XiX' invariant. Then 
the term A < >< x'l* > is a product of over- 
laps < if'ilxi >< x'jlVj > and the Hartree-Fock equa- 
tions read 

iVi - h)\ipi > = Xi\xi >< x!i\Vi > , (9a) 
<MVi-h) = X[<ip'i\xi><x'i\, (9b) 

Aj = reciprocal of the single-particle Green function. 



^^=K^{<X'^\{m-hr'\x^ >) 



(10) 



The operator h = t+U contains a self consistent mean 
field U, similar to the standard Hartree-Fock potential, 
and the propagation energies r]i are also given by self con- 
sistent formulae similar to those defining Hartree-Fock 
self energies. From (9), (10) we immediately read off the 
Hartree-Fock Hamiltonian corresponding to H', 



since there arc still infinitely many ways to adjust the 
phases and norms of Xi^Xi in such a way that X)X' 
main unchanged. We then extend the diagonalization 
of h' beyond the first N right and left cigcnstates, and 
obtain an infinite, biorthonormal set of orbitals ipa,ifi'i^, 
presumably complete. This set most likely includes con- 
tinuum states, to be normalized accordingly, hence 

< if'fflifa >= or S{r]0 - r]a), (14) 
{Va-h')\ip^>=0, < ip'f^\{r]p - h') = . (15) 

The equivalence of inversion and diagonalization on the 
many-particle level, based on the functionals i^(^', 5*) 
and i^'(5'',^), is reflected on the single-particle level: 
The homogeneous equations (9) for h' are strictly equiv- 
alent to the inhomogeneous TIMF equations, 

(% - h)\<fii >= \xi>, < ^'iliVi - h) =< x-l , (16) 

obtained from the functional F{^' with as 

Slater determinants. In particular the eigenvalues of 
/i'(A = 0) = /i, Eq.(15), coincide with the poles of D, 
where A = 0, when calculated from Eqs.(16). 

To improve the mean field approach by ph corre- 
lations, we introduce creation and annihilation opera- 
tors a\^,aa,ct = 1,2... 00 for orbitals fa', a^jifl,/? = 
1,2... 00 are the corresponding operators for ip^. Fi- 
nally a set of operators jj , Ji,i = 1,2 ... N is introduced 
for orbitals Xi, and a similar set J-\ for orbitals Xi- A 
second quantization representation of H' is then 

H' = Y.<'p'a\t\^fi> aia'fs 
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af3-yS 

-\-Ajlj2 • • ' J^^Jjs^ • ' • J2J1 • 



(17) 



In (17) one set only of operator pairs a]j,a^ appears. It 
anticommutes canonically according to (14). With | > 
as fermion vacuum we have 



iV 



N 



i$>=n«iio>, <$'i=<oin«i, 



(18) 



N 



AilXj ><x!i\ 



(11) 



bearing in mind (8). Standard manipulations of (9), to- 
gether with (8), show that <^i,V^ are biorthogonal. Fi- 
nally one can use normalizations such that 

< v'iWj >= h ; < f'i\xj >=< Xilvj >= lAj , (12) 

< X^\{V^ - hy^\X^ > 



V< XiliVi - h)-^\Xi > ' 



(13) 



as right and left quasi-particle vacua. Hence a complete 
set of ph operators is given by 



Bl 



at, a', B: 



m i' mi 



(19) 



for i = 1, 2 . . . A^, m = A^ -|- 1, . . . , 00, and the ansatz for 
boson operators on the RPA level is 

mi 

If one sets Y", = a priori, one falls back to the TDA. 
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While Wick's theorem is available to calculate the ma- 
trix elements of kinetic energy T and two-particle inter- 
action 1/, the matrix elements of A|x >< x'\ are less 
familiar. They are simplified, however, in the biorthogo- 
nal representation (12), and read, 

< "^'Wla'JX > - < ^'\X >< ^'JX^ > I < f'^\X^ >, 

(21a) 

< X'\ala',\^ > = < X'l* >< X-|V5m > / < X^\V'^ >, 
< >< ipWn\X^XJ > 



< 



^^\X^ >< ^'j\X] > 



(21b) 



/, t t ' '1^ <X'\'^>><X-X'j\^m^n> 

' <x'M ><x'j\vj > 

Let M be the number of ph pairs (mi) retained in a prac- 
tical calculation. The 2M solutions of the RPA equations 



(22) 



then split into 2 families of M solutions. For the first 
family, E* denotes the "excitation" energy which sepa- 
rates the RPA-correlated eigenstate > based on |$ > 
from another eigenstate of H' defined by the RPA ansatz 
>= Qiy\^ >■ In the second family, for which we 
will use a label V, —E* is the energy difference between 

< \iiL\ =<^'\Q- and the correlated state < *'| based on 

< $'1 . The solutions v of the first family can be told from 
those 17 of the second family by selecting as E* those RPA 
eigenvalues close to the eigenvalues of the TDA, which 
provides only M solutions. A straightforward, slightly 
tedious calculation, gives for (22) 

^ ' {C-nj.miX^^ + Ami,njY^^) = —E^Y^^ (23a) 
mi 



with 



Cnj.m'i. = - < iPiip'j\v\(prniPn > "A j 



< X'lX'jWni^Pn > 



Br, 



<x-m ><Xjm > 

(24) 

< ^'m^'JXiXj > 



A. 



< ^'M >< ^'j\Xj>' 
(25) 

_^j;^\<X^Mn><^^n]Xl>_ 
\ < X^\^^ >< V'jlXj > 

''^^ <x'M><^[\xi> /' ^ ^ 



using A = A<$'|x><x'|<l'> and the relation 

K=K=KI <^[\X^><X'^W^> ■ (27) 

Eqs.(23) show less symmetry than the usual RPA equa- 
tions 1^] , as expected from the non-hermitean H' as com- 
pared to H . However, B and C are symmetric under ex- 
change of the index pairs (mi) and (nj), the block matrix 



S = 



C A 
A B 



is symmetric and Eqs.(23) read 



S 



E*M ( J J , with M 



-1 

1 



(28) 



(29) 



It is then easy to show that the "symplectic diagonal- 
ization" problem (29) is equivalent to diagonalizing the 
antisymmetric matrix 



< $'| [[iJ',Q.],B^] I* > =E:< $'|[Q„B;^.]|$ > namely 



yv 



El-^S^I^ 



X" 
yv 



(30) 



(31) 



An elementary derivation of (31) from (29) involves the 
definition of any symmetric square root of 5* and its in- 
version, and the assumption that / 0. While trans- 
lational and rotational symmetry of H give rise to null 
eigenvalues in usual RPA, the \x >< x'\ -term in H' 
destroys these symmetries, hence eigenvalues E^ — are 
usually not to be expected. Special proofs are available 
when an exceptional incident occurs. 

The fact that A is antisymmetric proves that the spec- 
trum splits into two families of opposite eigenvalues, as 
announced. Consider now the AI eigenvalues Eq^,, com- 
ing from the diagonalization of A, and define the auxil- 
iary " perturbative RPA" matrix 



eC A 
A eB 



< e < 1. 



(32) 



with corresponding Ac- It is clear that, as e — > 0, half of 
the eigenvalues of the auxiliary problem (32) converge to- 
wards i?Q^, while the other M eigenvalues tend to ~Eq^. 
By continuity when e — > 1, this makes it easy to sepa- 
rate the eigenvalues of (29) into the two announced fami- 
lies, both being complex due to the nonhermiticity of H' . 
With the first family, the usual bosonic interpretation of 
the RPA-like Hamiltonian finally provides the correction 



Az(A) = 




(33) 



A function inversion of the final result Z2{h) = 2:9 (A) -|- 
Az(A) , generating a function K{z2), then provides the 
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improvement on the Green function, beyond the mean 
field, by inclusion of ph correlations. 

In conclusion, the key element of our study is the equiv- 
alence of the following two problems: 

1. ) Calculate the Green function D{z) at some fixed 
(complex) energy z by operator inversion of the exact 
Hamiltonian H. One solves an inhomogeneous problem 
with channel wave functions %, x' as inhomogeneity. 

2. ) Diagonalize a Hamiltonian H' which differs from H 
by the operator A|x >< comprising the reaction 
boundary conditions. In this homogeneous approach, one 
calculates the (complex) eigenvalues z as function of the 
parameter A = 1/D and obtains the desired Green func- 
tion by inversion of the function z{A). 

Thus operator inversion is replaced by operator diagonal- 
ization followed by function inversion. While problems 
1.) and 2.) can equally well be solved in the mean field 
approximation, we use 2.) to incorporate correlations in 
the wave function in terms of ph bosons. Although the re- 
sulting boson Hamiltonian is less symmetric than in stan- 
dard (hermitean) RPA, the remaining symmetry is suffi- 
cient to diagonalize the (non-hermitean) boson Hamilto- 
nian by canonical, but non-unitary transformation and 
to calculate z(A) in this boson approximation. 

We tested numerically the above method for a one- 
dimensional, separable, two-body Hamiltonian of relative 
motion 

< q'Q'lnlqQ >= (34) 
6{Q ~ Q') [qH{q - q') - \v\q' eM~^^q^) , 

with g, Q the relative and total momentum, which gives 
an obvious solution Q for the exact amplitude D{z) of 
(3). Our TIMF and RPA numerics, however, used the 
single-particle representation which becomes mandatory 
for antisymmetrized, many-particle calculations. The 
channel waves x, x' were taken as boosted Gaussian func- 
tions. Mass and potential parameters as well as boost 
and width of x, %' were chosen to be typical for nuclear 
physics problems; for technical details we refer to ||]. The 
RPA approximation of D{z), based on (23) and (33), in- 
deed improved the TIMF results (see Fig.l) obtained by 
solving (16). 
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FIG. 1. Real and imaginary part of the exact, TIMF- 
and RPA-amplitudes vs. Re E, at Im E=2.0 for a solvable 
two-body model with gaussian interaction. 
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